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Abstract 



An infinite system of neutrons interacting by a model pair potential is considered. 
We investigate a case when this potential is sufficiently strong attractive, so that its 
scattering length a tends to infinity, a — > — oo. It appeared, that if the structure of 
the potential is simple enough, including no finite parameters, reliable evidences can 
, be presented that such a system is completely unstable at any finite density. The 

I incompressibility as a function of the density is negative, reaching zero value when the 

' density tends to zero. If the potential contains a sufficiently strong repulsive core then 

. the system possesses an equilibrium density. The main features of a theory describing 

such systems are considered. 



^ , PACS. 24.10.Cn Many-body theory - 21.10.Dr Binding energies 



There exists a well known problem in many body physics related to the description of the 



■ ground state properties of an infinite system composed of interacting fermions. In general, 

this description is based usually on tedious numerical calculations, particularly when the 
interaction is rather strong. The well known exceptions from this situation, when it is possible 
to calculate the ground state properties analytically, are the Random Phase Approximation 
(RPA) for a high density electron gas ^ and the low density approximation for dilute gases 
0. In both cases the kinetic energy Tk is much bigger then the interaction energy Eint 
of the system. This allows to apply some kind of a perturbation theory. In the case of 
homogeneous electron liquid it turns out that the analytical RPA-like description is also 
possible not only at very high but medium densities when Tk ~ Eint H, 3 • Similar extension 
of the range of validity is impossible in the case of fermion systems at low densities p: there 
the gas approximation is not applicable if Tk ~ Eint- In the cases when the pair interaction 
is attractive and sufficiently strong, the system can have a quasi equilibrium or equilibrium 
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states in which ~ -Em*- On the other hand, these states are preceded by special points 
with density p values at which the incompressibility K{p) of the system tends to zero. Thus, 
if it would be possible to predict the existence of such points then in principle it would 
become possible to conclude that the system has at least a quasi equilibrium state. 

In this Short note we address the above mentioned problem and consider the ground 
state properties of the infinitely extended multi-fermion system. We demonstrate that it can 
be done analytically provided that the pair interaction between fermions is characterized 
only by the scattering length a — oo. One can say in this case that the scattering length 
is the dominant parameter of the problem under consideration. Such an investigation is 
of great importance since it can be applied to fermion systems interacting via potentials 
with not only infinite, but also sufficiently large a. For instance, the scattering length 
a of neutron-neutron interaction is about —20 fm, thus greatly surpassing the radius of 
the interaction tq. On the other hand, it is possible now to prepare artificially a system 
composed of Fermi atoms interacting by an artificially constructed potential with almost 
any desirable scattering length, similarly to that how it is done for the trapped Bose gases, 
see e.g. An experimental study, performed on such Fermi-system would be of great 

importance presenting new information on the behavior of dilute gases and the gas-liquid 
phase transition. 

Let us consider the interaction of two isolated particles. We assume that this interaction is 
of finite radius tq, which is small, so that ppf'o 1 (pf = (Svr^p)^/^ is the Fermi momentum), 
but its strength is such that the scattering length is negative and infinitely big, a — oo. 
We assume also, that the density p of the system under consideration is homogeneous. As 
it will be demonstrated below, in such a case the system is located in the vicinity of a 
phase transition, transforming it into a strongly correlated one. Therefore, the problem of 
calculating its ground state properties has to be treated for the most part qualitatively. 

Let us start considering general properties of a Fermi system with some attractive two- 
particle bare interaction V{r), which is sufficiently weak to create a two-particle bound state. 
Assume, that the scattering length a corresponding to this potential is negative and finite. 
The ground state energy density E{p) can be approximated by Skyrme-like expression 0, 

E{p) = ^P + t,p^ + hpy\ (1) 

Here M is the particle mass, and p is their density in the system. The first term of eq. 
(1) is the kinetic energy T^, while the second and the third are related to the interaction 
energy Eint determined by the potential V{r). The second term which is proportional to to 
gives a proper description in the gas limit. The third term provides the behavior of E{p) at 
higher densities, including that of the equilibrium density. Such structure of E{p) appears 
if the interaction is sufficiently attractive so that to < 0, and > 0. Note, that eq. (1) 



2 



presents at least a qualitative description of the system under consideration giving a rather 
simple and reasonable presentation of the function E{p). A more precise picture of the 
energy dependence upon density can be obtained using more sophisticated functionals for 
the ground state energy |^. 

If the potential V is of short range and purely attractive, then in the Hartree-Fock 
approximation the ground state energy Ehf is given by the following expression 

3„2 

EHFip) = J^P + tHFP^ (2) 

where the parameter to = tHF , being negative, is entirely determined by the potential 
V{r). For instance, in the case of a short range 5-type interaction one has Ihf = —vo/^, 
with vq being the corresponding strength of the potential. Eq. (2) shows that at small 
densities Ehf > due to the kinetic energy term, but at sufficiently high densities p ^ oo 
the Hartree-Fock energy becomes dominating, leading to the collapse of the system, with 
Ehf — oo. Keeping in mind that the Hartree-Fock approximation gives the upper limit 
to the binding energy Ehf > E, one can conclude that the system does not have, in this 
case, an equilibrium density pe and energy Ee since Eg — >• — oo when p — > cxd 0. Note, that 
for a given and finite total number of particles N, the HF energy is not going to infinity and 
the system collapses into a small volume with the radius tq, with the density p ~ N/vq. It 
is evident that the function E{p) is positive at small densities, if the parameter to is finite. 
Therefore, it must have at least one maximum at the density pm before it becomes negative, 
on the way to E ^ — oo. If the potential V{r) includes a repulsive core at sufficiently short 
distances, then > . As a result, the system has an equilibrium density and energy, pe 
and Ee, respectively, determined by the repulsive core strength and its radius Vc ~ tq. 

Now let us apply eq. (1) to demonstrate the most important features of the system under 
consideration: 

a) when p — the third term on r.h.s. in eq. (1) can be omitted. The kinetic energy is 
relatively very big, 3> Eint, and to ~ cl, with a < being the scattering length. In that 
case we have a dilute Fermi gas with positive pressure P and incompressibility K, the latter 
being determined by the equation, see e.g. 0, 

^ dE^p) 

K{p) = p (3) 

b) on the way to higher densities, which can be achieved by applying an external pressure, 
the system reaches the density pci < Pm at which the incompressibility is equal to zero, 
K{pci) = 0. Remembering that at the maximum the second derivative is negative, one can 
conclude, as it is seen from eq. (3), that K{pm) < 0. In the range Pc2 > P > Pci the 
incompressibility is negative, i^' < 0, and as a result the system becomes totally unstable. In 
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fact, in this density range all calculations of the ground state energy are meaningless since 
such a system cannot exist and thus be observed experimentally | T0[] ; 



c) at some point p = Pc2 > Pm the contribution due to the repulsive core becomes 
sufficiently strong to prevent the further collapse of the system. The incompressibility attains 
K = at Pc2 < Pe, being positive at the higher densities. Finally, the system becomes stable 
at p > Pc2, reaching equilibrium density at pe with equilibrium energy equal to i?e- It is 
obvious that K{pe) > being proportional to the second derivative at the minimum, see eq. 
(3). It should be kept in mind that in this density domain, p > pc2, the function E{p) is 
determined by the repulsive part of the potential which makes ts > 0. As it was mentioned 
above, without this component of V{r) the system's energy would infinitely increase, Ehf 
— oo, with density growth, p — oo, thus inevitably collapsing. 

One could expect in principle the existence of metastable states at p > pd if the potential 
V{r), even being pure attractive, would have a complicated structure. It can be said that 
there could exist parameters of V{r), which are able to open the possibility for the metastable 
states to be formed. On the other hand, a system of fermions interacting via a short- 
range, finite scattering-length, 5-type potential Vs , can be stable only in the dilute gas 
regime. While at the densities p > pd the incompressibility K becomes negative, the system 
collapses. Indeed, the potential Vs has no structure to ensure any metastable states at the 
densities p > Pci- As a result, one can write down a dimensionless expression for the ground 
state energy as a function of the only variable z = ppa 0, 



aEiz) = z'il + Piz)), (4) 
with a = lOvr^Ma^. In the low density limit, \apF\ <^ 1 and when the interaction has the 



radius ro, eq. (4) reads |jTT 



aE{z) = z' 



97r ilTT^ V a / V a 



(5) 



Here the function 'j{y,z) is of the order of one, 7(1/, 2;) ~ 1. It is seen from eq. (5) that as 
soon as the scattering length becomes big enough, \a\ ^ tq, one can omit the contribution 
coming from the function 7 and neglect all the term proportional to (rg/a)^. Then eq. (5) 
reduces to eq. (4). Thus, in the case \a\ — > cx) we can use eq. (4) to determine the ground 
state energy E. Eq. (4) is valid up to the density pd which is a singular point of the function 
(3{z), since beyond this point K < 0, and the system is completely unstable. On the other 
hand, there is no physical reasons to have another irregular point in the region < p < pd. 
We note, that eq. (5) is as well valid up to its own density p^^^ ^ pd, provided \a\ ^ tq. 
Using eq. (3) for the incompressibility, one can calculate the position of the point Zd where 
K = 0. Denoting the corresponding z as Zd = cq, where cq is a dimensionless number. 
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one is led to the conclusion that pci 



provided a is sufficiently large to be the only 



dominating parameter. The system has only one stable region at small densities p < pci 
which decreases and even vanishes as soon as a — > — oo. One could expect that |co| oo 
as soon as a becomes the dominant parameter so that the above given expression for f3{z) is 
valid in the whole domain \z\ < oo. On the other hand, there exists another singular point 

in the function E{p) and the position of this point which corresponds to Pc2 depends 
mainly on the parameters such as tq, core radius Vc of V{r) rather than on the scattering 
length a. As to the function f3{z), by definition it does not contain any information about 
Pc2- Therefore, in order to take into account the existence of pc2, one has to recognize that cq 
is finite, and the densities pci and pc2 are different. As a result, the function (3 is determined 
in fact only in the region \z\ < \zci\. 

Now let us consider the calculation of the ground state energy E of the system when the 
density approaches p — > Pci ~ |a. 

Pc2 are missed in calculations because of the lack of the self consistency 
relates the linear response function of system with its incompressibility K 



^ from the low density side. As a rule, points pd and 



which 



Xiq 0,iiJ ^ 0) = - 



(6) 



As we shall see below, these points can give important contributions to the ground state 
energy. To see it we express the energy of a system in the following form (see e.g. 



E{p) 



(7) 



where x(g, ict;, (?) is the linear response function on the imaginary axis and v{q) = gV{q), 
with V{q) being the Fourier image of V{r). The integration over u goes from — oo to +oo, 
while the integration over the coupling constant g runs from zero to the real value of the 
coupling constant, i.e. to g = 1. At the point p = pd the linear response function has a 
pole at the origin of coordinates q = 0, u = due to eq. (6). At the densities p > pd 
the function iuj) has poles at finite values of the momentum q and frequencies iu. This 
prevents the integration over iuj, making the integral in eq. (7) divergent. Thus, we conclude 
that it is the contribution of these poles that reflects the system's instability in the density 
range Pd < p < Pc2- Note, that violations of eq. (6) lead to serious errors in the calculation 
of the ground state energy. Eq. (7) can be rewritten, explicitly accounting for the effective 
interparticle interaction R{q,iu!,g), (see e.g. W^), in the following form 



E{p) =Tk + EH-\j 



+ 27ipS{uj) 



v{q) 



dq duj dg 
7(2^- 
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Here xo is the linear response function of noninteracting particles, while x is given by the 
following equation 

x{q,^) = -. — ^7 — ^ — 7 — (9) 

1 - R{q,uj)xo{q,uj) 

It is seen from eqs. (6) and (9) that the denominator {1 — Rxo) vanishes at p ^ pd while the 



radius of correlation tends to infinity |]I0[. Thus, it is impossible to present the denominator 
as a power series in Rxo approximating the expansion by the finite number of terms. This 
result is quite obvious since pd is a singular point in the function E{p) which makes it 
impossible to expand that function in the vicinity of this point. Therefore, one should try to 
satisfy eq. (6) in order to get proper results for the ground state calculations in the vicinity 
of the instability points. Such an approach was suggested in ^, [T^ and is based on the 
exact functional equation for the effective interaction R{q,u!,g), 

2 6p\q,uj)J 1 - R{k,iw,g}xo{k,tw} g{2Tcp 

As a result, the linear response function x given by eq. (9) automatically satisfies eq. (6) 
I^ . Our preliminary calculations |l5l, based on eq. (9) and applied to the case when 



the scattering length is sufficiently large but finite, confirm the result that pci ~ kl~ • 

Let us suppose for a while that the bare potential is pure attractive. Then, the interval 
of the densities [0,pci] within which the system is stable vanishes with the growth of \a\. As 
a result, in the limit a = —oo the incompressibility becomes negative < 0, making the 
considered system completely unstable at any density. Thus, the point at which a = — oo 
is the only point of the system's instability at all the densities. As soon as the scattering 
length deviates from its infinite value, that is +oo > a > — oo the system comes back to 
its stable state at list in the range of the density values p < pd ~ |a|~^. It is of interest to 
understand whether it is possible to prove by e.g. numerical calculations, that pd ~ 
when a — oo. From our point of view, at least at this moment, the answer is "no". We are 
dealing with a system located in the vicinity of a phase transition, which transforms it into 
a strongly correlated one. As a result, it is hard to believe that the numerical calculations 
could be reliable. On the other hand, it is not really necessary to carry out numerical 
calculations if the problem allows a qualitative analysis. It was argued above, that there 
exists the only parameter to characterize the system which is the scattering length a. In fact 
the scattering length determines only the specific point pd at which the incompressibility 
vanishes, separating the region of a dilute gas from the region of the system's instability. As 
soon as a — » — oo this last and the only parameters vanishes, driving the point pd ~ lo-l^^ 
of the curve E{p) to the origin of coordinates. Thereafter, the system becomes unstable at 
all densities. And vice versa, as soon as the scattering length becomes finite the system is 
stable at list within the interval p < Pd ~ kl"^- 
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Note, that as it follows from our consideration, any Fermi system possesses an equilibrium 
density and energy if the bare particle-particle interaction contains a repulsive core and its 
attractive part is strong enough, so that a —oo. Indeed, at sufficiently small densities 
the ground state energy is negative (since the incompressibility K < 0) and the system will 
collapse until the core stops the density growth. Therefore, the minimal value of the ground 
state energy must be negative when the repulsive core will enter the play to prevent the 
system from the further collapse. It is worth to remark, that superfluid correlations cannot 
stop the system squeezing, since their contribution to the ground state energy being negative 
increases in the absolute value with the growth of the density. 

A liquid similar to the model one considered in this paper exists in Nature. This is liquid 
^He. If a helium dimer exists, its bound energy does not exceed 10~^ meV while the ground 



state energy of hehum liquid is about 2 * 10~ meV per atom Because of this huge 

difference in binding energies, it is evident that there is no essential contribution coming 
from the binding energy of the dimer to the ground state energy of the liquid. In fact, the 
numerical calculations show that the pair potential is rather weak to produce the dimer He2 
T6[| . Thus, one can reliably consider an infinite homogeneous system of Helium atoms as 



consisting of particles interacting via pair potential, characterized by a very big but finite 
scattering length \a\ ^ tq. Let us make also the following additional remark. It seems 
quite probable that the neutron-neutron scattering length (a ~ —20 fm) is sufficiently large 
to permit the neutron matter to have an equilibrium energy and density fl^. Therefore, 
calculations of a neutron matter satisfying eq. (7) are quite desirable. 

In summary, the homogeneous system of interacting fermions was considered. It was 
shown that when the scattering length a is negative and sufficiently large the fermion matter 
becomes a strongly correlated system at the densities p ~ \a\~^. Therefore, the consideration 
of such a system is connected to a number of problems which yet persist and have to be 
resolved. At the same time, the qualitative consideration presented above gives strong 
evidences that the point pd at which the incompressibility vanishes is defined by pci ~ 
\a\~^ provided the scattering length is the dominant parameter of the problem. Thus, a 
homogeneous system composed of fermions, interacting via a pure attractive potential, at 
a ^ — oo is completely unstable at all the densities, with the incompressibility as a function of 
the density being always negative. As soon as the density p goes to zero the incompressibility 
goes to zero as well. 

We thank G.F. Bertsch for attracting our attention to the discussed above many body 
problem which he has in fact suggested One of us (VRS) is grateful to A. Bulgac 

for valuable discussions, and to the Department of Physics of the University of Washington 
where part of this work was done, for hospitality. This research was funded in part by INTAS 
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